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Summary. The evolution of the probability distributions of Japan and US major 
market indices, NIKKEI 225 and NASDAQ composite index, and JPY/DEM and 
DEM /U SD currency exchange rates is described by means of the Fokker-Planck 
equation (FPE) . In order to distinguish and quantify the deterministic and random 
influences on these financial time series we perform a statistical analysis of their 
increments Ax(A(t)) distribution functions for different time lags A(t). From the 
probability distribution functions at various A(t), the Fokker-Planck equation for 
p(Ax(t), A(t)) is explicitly derived. It is written in terms of a drift and a diffusion 
coefficient. The Kramers-Moyal coefficients, are estimated and found to have a simple 
analytical form, thus leading to a simple physical interpretation for both drift 
and diffusion coefficients. The Markov nature of the indices and exchange rates 
is shown and an apparent difference in the NASDAQ is pointed out. 
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1 Introduction 

Recent studies have shown that the power spectrum of the stock market fluc- 
tuations is inversely proportional to the frequency on some power, which 
points to self-similarity in time for processes underlying the market [1, 2]. 
Our knowledge of the random and/or deterministic character of those pro- 
cesses is however limited. One rigorous way to sort out the noise from the 
deterministic components is to examine in details correlations at different 
scales through the so called master equation, i.e. the Fokker-Planck equation 
(and the subsequent Langevin equation) for the probability distribution func- 
tion (pdf) of signal increments [3] . This theoretical approach, so called solving 
the inverse problem, based on rigorous statistical principles [4, 5], is often 
the first step in sorting out the best model(s). In this paper we derive FPE, 
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directly from the experimental data of two financial indices and two exchange 
rates series, in terms of a drift and a diffusion coefficient. We would 
like to emphasize that the method is model independent. The technique al- 
lows examination of long and short time scales on the same footing. The so 
found analytical form of both drift and diffusion D^ coefficients has 
a simple physical interpretation, reflecting the influence of the deterministic 
and random forces on the examined market dynamics processes. Placed into 
a Langevin equation, they could allow for some first step forecasting. 

2 Data 

We consider the daily closing price x(t) of two major financial indices, 
NIKKEI 225 for Japan and NASDAQ composite for US, and daily ex- 
change rates involving currencies of Japan, US and Europe, JPY/DEM 
and DEM /USD from January 1, 1985 to May 31, 2002. Data series of 
NIKKEI 225 (4282 data points) and NASDAQ composite (4395 data points) 
and are downloaded from the Yahoo web site (http : / / finance.yahoo.com/). 
The exchange rates of JPY/DEM and DEM /USD are downloaded from 
http : / /pacific.commerce.ubc.ca/xr/ and both consists of 4401 data points 
each. Data are plotted in Fig. l(a-d). The DEM /USD case was studied in 
[3] for the 1992-1993 years. See also [6], [8-10] and [11] for some related work 
and results on such time series signals, some on high frequency data, and for 
different time spans. 



3 Results and discussion 

To examine the fluctuations of the time series at different time delays (or 
time lags) At we study the distribution of the increments Ax = x(t + At) — 
x(t). Therefore, we can analyze the fluctuations at long and short time scales 
on the same footing. Results for the probability distribution functions (pdf) 
p(Ax, At) are plotted in Fig. 2(a-d). Note that while the pdf of one day time 
delays (circles) for all time series studied have similar shapes, the pdf for 
longer time delays shows fat tails as in [2] of the same type for NIKKEI 225, 
JPY/DEM and DEM /USD, but is different from the pdf for NASDAQ. 

More information about the correlations present in the time series is given 
by joint pdf's, that depend on N variables, i.e. p N (Ax±, At\\ Axn, Atpj). 
We started to address this issue by determining the properties of the joint 
pdf for N = 2, i.e. p(Ax 2 , At 2 ; Ax\, At\). The symmetrically tilted character 
of the joint pdf contour levels (Fig. 3(a-c)) around an inertia axis with slope 
1/2 points out to the statistical dependence, i.e. a correlation, between the 
increments in all examined time series. 

The conditional probability function is 
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Fig. 1. Daily closing price of (a) NIKKEI 225, (b) NASDAQ, (c) JPY/DEM and 
(d) DEM/USD exchange rates for the period from Jan. 01, 1985 till May 31, 2002 



p (Ax i+ i,At i+ i \Axi,Ati ) = — — (1 

p(Axi,Ati) 

for i — 1, N — 1. For any At 2 < Ati < Ati, the Chapman-Kolmogorov 
equation is a necessary condition of a Markov process, one without memory 
but governed by probabilistic conditions 



p(Ax 2 ,At 2 \Axi,Ati) = J d(Axi)p(Ax2,At2\Axi,Ati)p(Axi,At i \Ax 1 ,Ati). 

(2) 

The Chapman-Kolmogorov equation when formulated in dif ferential 
form yields a master equation, which can take the form of a Fokker-Planck 
equation [4]. For r = log 2 (32/At), 
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in terms of a drift D^(Ax,t) and a diffusion coefficient D^(Ax,t) (thus 
values of r represent Ati, i = 1, ■■■)• 
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Fig. 2. Probability distribution function p(Ax,At) of (a) NIKKEI 225, (b) NAS- 
DAQ, (c) JPY/DEM and (d) DEM/USD from Jan. 01, 1985 till May 31, 2002 
for different delay times. Each pdf is displaced vertically to enhance the tail behav- 
ior; symbols and the time lags At are in the insets. The discretisation step of the 
histogram is (a) 200, (b) 27, (c) 0.1 and (d) 0.008 respectively 



The coefficient functional dependence can be estimated directly from the 
moments M^ k ' (known as Kramers-Moyal coefficients) of the conditional prob- 
ability distributions: 

M (k) = — [ dAx'{Ax - Ax) k p(Ax ,t + At\Ax,t) (4) 
At J 

dW(Ax,t) = ^limAfW (5) 
re! 

for At — > 0. The functional dependence of the drift and diffusion coefficients 
and for the normalized increments Ax is well represented by a line 
and a parabola, respectively. The values of the polynomial coefficients are 
summarized in Table 1 and Fig. 4. 

The leading coefficient (oi) of the linear dependence has approxi- 

mately the same values for all studied signals, thus the same deterministic 
noise (drift coefficient). Note that the leading term (62) of the functional de- 
pendence of diffusion coefficient of the NASDAQ closing price signal is about 
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Fig. 3. Typical contour plots of the joint probability density function 
p(Ax 2 ,At 2 ;Ax 1 ,At 1 ) of (a) NIKKEI 225, (b) NASDAQ closing price signal and 
(c) JPY/DEM and (d) DEM /USD exchange rates for At 2 = 1 day and 
Ati = 3 days. Contour levels correspond to logioP = —1.5, —2.0, —2.5, —3.0, —3.5 
from center to border 



Table 1. Values of the polynomial coefficients defining the linear and quadratic 
dependence of the drift and diffusion coefficients D' 1 ' = ai(Ax/a) + ao and D^ = 
b2(Ax/a) 2 +bi(Ax/a)+bo for the FPE (3) of the normalized data series; a represents 
the normalization constant equal to the standard deviation of the At— 32 days pdf 





ffli 


ao 


b 2 


bi 


bo 


a 


NIKKEI 225 


-0.54 


0.002 


0.18 


-0.004 


0.003 


1557.0 


NASDAQ 


-0.49 


-0.0004 


0.30 


-0.010 


0.001 


198.11 


JPY/DEM 


-0.55 


0.002 


0.17 


-0.002 


0.004 


2.9111 


DEM/USD 


-0.49 


-0.001 


0.16 


-0.004 


0.004 


0.0808 



twice the leading, i.e. second order coefficient, of the other three series of 
interest. This can be interpreted as if the stochastic component (diffusion 
coefficient) of the dynamics of NASDAQ is twice larger than the stochastic 
components of NIKKEI 225, JPY/DEM and DEM /USD. A possible rea- 
son for such a behavior may be related to the transaction procedure on the 
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Fig. 4. Functional dependence of the drift and diffusion coefficients D (1 ^ and 
for the pdf evolution equation (3); Ax is normalized with respect to the value of the 
standard deviation a of the pdf increments at delay time 32 days: (a,b) NIKKEI 225 
and (c,d) NASDAQ closing price signal, (e,f) JPY/DEM and (g,h) DEM /USD 
exchange rates 
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Fig. 5. Equal probability contour plots of the conditional pdf p(Ax2, At2\Axi, Ati) 
for two values of At, At\ = 8 days, AI2 = 1 day for NASDAQ. Contour levels 
correspond to logiop=-0. 5,-1. 0,-1. 5,-2. 0,-2. 5 from center to border; data (solid line) 
and solution of the Chapman Kolmogorov equation integration (dotted line); (b) 
and (c) data (circles) and solution of the Chapman Kolmogorov equation integration 
(plusses) for the corresponding pdf at Ax2 = -50 and +50 



NASDAQ. Our numerical result agrees with that of ref. [3] if a factor of ten 
is corrected in the latter ref. for 62 • 

The validity of the Chapman-Kolmogorov equation has also been verified. 
A comparison of the directly evaluated conditional pdf with the numerical 
integration result (2) indicates that both pdf 's are statistically identical. The 
more pronounced peak for the NASDAQ is recovered (see Fig. 5). An analyt- 
ical form for the pdf's has been obtained by other authors [6, 10] but with 
models different from more classical ones [8]. 



4 Conclusion 

The present study of the evolution of Japan and US stock as well as foreign 
currency exchange markets has allowed us to point out the existence of deter- 
ministic and stochastic influences. Our results confirm those for high frequency 
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(1 year long) data [7, 10]. The Markovian nature of the process governing the 
pdf evolution is confirmed for such long range data as in [3, 7, 8] for high 
frequency data. We found that the stochastic component (expressed through 
the diffusion coefficient) for NASDAQ is substantially larger (twice) than for 
NIKKEI 225, JPY/DEM and DEM/USD. This could be attributed to the 
electronic nature of executing transactions on NASDAQ, therefore to different 
stochastic forces for the market dynamics. 
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